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Qh We present CP asymmetry measurements using the full BABAR dataset 

of 471 million BB pairs for the charmless B decays: B° — > p°K*°, B° — > 
f (980)K*°, B° -»■ p~K*+, and for the Dalitz plots of B + -»■ K + K~K + , 
B + K° S K° S K + and 5° -> K S K + R-. 
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1 Introduction 

in 

, J Charmless B decays are useful probes of the dynamics of weak and strong interactions. 

The interference between tree level and penguin contributions to the same final state 
can give rise to direct CP violation. The relative weak phase between tree and penguin 
amplitudes probes the Unitarity Triangle angle 7. In addition to this, enhancement 
in CP asymmetries with respect to the Standard Model expectations can signify 
beyond-the-Standard-Model particle contributions at loop level. 

In this paper, we present the results of analyses conducted by the BABAR collab- 
oration for the following charmless decay modes: B° — > p°K*°, B° — > fo(980)K*°, 
B° -> p~K*+, B+ -> K+K-K+, B+ -> K° S K° S K+ and B° -> K S K + K~. These 
analyses make use of the full BABAR Y(4S) dataset. This consists of 467 million 
BB pairs collected by the BABAR detector at the PEP-II B factory, which collides 
e + e~ asymmetric-energy beams at the T(AS) resonance pQ. The B meson candidates 
are characterised by using two kinematic variables. We take advantage of the pre- 



cise kinematic information from the beams to form the variables mgs = \/f — Vb 
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and AE = E* B — where {Eg,pg) is the B meson four-momentum in the e + e~ 
centre-of-mass frame and \fs is the centre-of-mass energy. Signal events are expected 
to peak around the B mass for mgs and around zero for AE. There are two main 
types of backgrounds: the very abundant qq background, where q is either a u, , s 
or ^quark, and the background arising from B decays to other final states. To distin- 
guish B meson candidates from the continuum qq background, variables describing 
the topology of the event are combined in a multivariate analyser (MVA), such as a 
neural network or a Fisher discriminant, in order to maximise their discriminating 
power. The variables mss, AE and the output of the MVA can either have selection 
requirements placed upon them or be supplied as inputs to a maximum likelihood 
fit. BB backgrounds are reduced by vetoing events containing a candidate consistent 
with a charm particle, and identifying remaining events by adding a category in the 
maximum-likelihood fit. 

2 Direct CP violation in B° p°K*°, B° / (980)1T ( 
and B° -> p~K* + 

The measurement of direct CP violation in B° — > pK* can be used to extract the CKM 
angle 7 from b — » s penguin contributions [2] . However CP asymmetries in these decay 
modes have not been previously well measured due to limited statistical sensitivity and 
predictions for these decays are not robust. BABAR published measurements of direct 
CP violation in B° ->■ p°K*° (892) and B° ->■ f K*° (892) based on approximately 
half the BABAR dataset. CP asymmetries were found to be consistent with no direct 
CP violation. However this analysis revealed only 3.5cr evidence for the decay B° — > 
f K*° (892), and found no statistically significant rate for B° p~K*+ (892) [3]. 

Neutral B candidates are reconstructed from their decays to p° or / candidate 
decaying to 7r + 7r~, or as a p~ candidate decaying to 7r~7r°. These are combined with a 
K*° or K* + candidates reconstructed as decays to K + tt~ or K + n°, respectively. The 
signal region is defined in terms of rriEs, AE and a range in the Kir invariant mass 
around the K*(892) mass. To account for resolution effects due to the 7r° meson, 
the AE signal region for the p~K* + decay is defined to be asymmetric around zero, 
namely —0.17 < AE < 0.1 GeV. In addition to a Fisher discriminant formed from 
four event-shape variables, continuum background is also suppressed by applying 
a selection on the angle between the thrust axis of the B candidate and that of 
the rest of the event. A main source of B backgrounds for the p°/foK*° decays 
comes from B° — >■ f2(l270)K*° . To identify the /2(1270) contribution in data, a 
maximum-likelihood fit using m-ES; AE, the Fisher discriminant and Kir invariant 
mass is performed. s Weights [1] are used to obtain the tttt invariant mass for the 
K*° signal events. A fit to the distribution gives a total of 47 ± 3 /2(1270) events in 
the signal region. 
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The final maximum-likelihood fit includes seven observables: /bes, AE, the Fisher 
discriminant, Kir and it + it~ invariant masses, and two helicity angles. Figure [l] shows 
the results of the fit to data. This is the first observation of B° — > foK*° (892) and 
B° — > p~ K* + (892) with significance greater than 5a. Results for CP asymmetries are 
listed in Table [T} Systematic uncertainties on CP asymmetries include a potential 
bias due to the difference in the interactions of the K + and K~ with the detector 
material (~ 1%). No significant direct CP violation is observed in any of the three 
decay modes [3]. 
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Figure 1: m^s projections (a,c) and tt + tt invariant mass projections (b,d) for the 
p°/f K*° (892) fit (a,b) and the p~K* + (892) fit (c,d). The points with error bars show 
the data distribution, the solid curve shows the overall fit result, the dot-dashed curve 
shows total background, and the dashed curve the total signal contribution. For the 
K*° mode, the overall signal has been split into a dashed p° component and a dotted 
f component. 



Table 1: The measured CP asymmetries for the p°/foK*° and p K* + decay modes. 
Decay mode A C p 
B u ->• p°K* {) (892) -0.06 ±0.09 ±0.02 

B° -»■ / (980)ir*°(892) 0.07 ± 0.10 ± 0.02 
B° -> p-K* + (892) 0.21 ±0.15 ±0.02 



3 Dalitz plot analyses of B + K + K~K + , B + -> 
K° S K° S K + and B° -> K S K + K~ 

The decays B + -»■ K + K~K + , B+ -»■ K° S K° S K + and B° K° S K + K~ are all b -»■ s 
penguin-dominated decay modes. The inclusive direct CP asymmetries for these de- 
cays are predicted to be small, and in the range (0.0 - 4.7)% [HI [7j. Any enhancement 
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of the CP asymmetry values would be an indicator of new physics. Previous analyses 
of these decay modes have revealed complex Dalitz plot structures that include a very 
broad nonresonant contribution, and a poorly understood scalar resonance, dubbed 
the /x(1500) [SI [HI UUl HTj . An example Dalitz plot Monte-Carlo simulation model for 
B + K + K~K + , which includes an /x(1500) scalar resonance, is shown in Figure [2] 
To create a correct Dalitz plot model, it is necessary to understand these structures. 
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Figure 2: Dalitz plot model from K + K K + MC including an /x(1500) component. 

Events are reconstructed as three charged-particle tracks consistent with the kaon 
hypothesis for B + — > K + K~ K + , two charged-particle tracks and a K® — > tt + tt~ or 
7r°7r° candidate for B° — > K®K + K~ , or one charged-particle track and two K® — > 
7r + 7r~ candidates for B + — > K° S K®K + . This analysis uses Legendre polynomial mo- 
ments to test the Dalitz plot model of each B — > KKK decay mode. Two different 
models are studied: model A includes a spin-zero /x(1500) component and model B 
replaces the /x(1500) with / (1500) and f' 2 {1^2^) components. The data are fitted 
by varying the mass and width of the /x(1500) and using an exponential S-wave term 
for the nonresonant component for model A, and a polynomial nonresonant model 
with S-wave and P-wave terms for model B. The Legendre polynomial moments are 
defined in terms of KK mass intervals as: 

N 

(Pi(cos6 3 ,m KK )) &J2 P i( cosd i)> C 1 ) 
i=i 

where 83 is the KK helicity angle, and iV is the number of events in the mxK mass 
interval. The s Weights for signal moments are then plotted as a function of KK 
invariant mass, one example of which is shown in Figure [3j with the distributions for 
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model A and B overlaid. For all three decay modes, it was found that a model includ- 
ing the /x(1500) does not describe the data distribution well. However a significant 
improvement in log likelihood is achieved when using model B. 
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Figure 3: The B + — > K + K~K + second-order angular moments for signal- weighted 
data (gWeights) in the region m K + K -i ow > 1.04GeV/c 2 compared to model A 
(dashed) and model B (solid). 

Tableggives the CP asymmetry results for B + -»■ K + R-R + and B° K° S K + K~. 
In the three-charged-kaon decay mode, the CP asymmetry for 0(1O2O).K' + differs from 
zero by about 2.8 standard deviations. Experimental precision for the CP asymme- 
tries in B° — > K Q S K + K~ is however not as good as for B + — > K + K~ K + . Log 
likelihood plots for the CP asymmetries for 0(1O2O)K + and 0(1020)^° are shown in 
Figure [4j The maximum- likelihood fit for B + — > K®K®K + gives 15 different solutions 
all within 3a of the global minimum. The amplitudes for resonances in this Dalitz 
plot are poorly constrained by the current dataset and therefore only an overall direct 
CP asymmetry is quoted, which is consistent with no direct CP violation [12] . 

Table 2: Results for the CP asymmetries for the resonances observed in the Dalitz 
plots for B+ -> K + K~K + and B° -> K° S K + R-. 



Decay mode 


A C p 


B+ -»■ 0(lO2O)fT + 
B+ f (980)K+ 
B + -»■ ^(1525)iT + 
B + ->■ K+K~K + (NR) 


0.13 ±0.04 ±0.01 
-0.08 ±0.08 ±0.04 
0.14 ±0.10 ±0.04 
0.06 ± 0.04 ± 0.02 


B° -»■ 0(1O2O)K° 
B° f (980)K° s 
Other 


-0.05 ±0.18 ±0.05 
-0.28 ±0.24 ±0.09 
-0.02 ±0.09 ±0.03 
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Figure 4: Scan of the 2A In C with (solid curve) and without (dashed curve) sys- 
tematic uncertainties, as a function of the CP asymmetry in 0(lO2O)fT + (left) and 
0(1020)^° (right). 



4 Conclusion 

BABAR is still actively producing new physics results for charmless hadronic B decay 
modes. Most of the direct CP violation results presented here are consistent with 
the Standard Model expectations. However BABAR measures the CP asymmetry 
of B + — > 0(lO2O)i^ + to be higher than the Standard Model expectation by 2.8cr. 
This is not yet evidence of new physics, and more data are required to confirm this 
measurement. Both LHCb and Belle 2 will be able to study the Dalitz plot structures 
of B + — > K + K~K + and B° — > K®K + K~ with much larger data samples. 
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